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ABSTRACT 

In recent years, accurate observational constraints become available for an increasing number of Galactic 
X-ray binaries. Together with proper motion measurements, we could reconstruct the full evolutionary history 
of X-ray binaries back to the time of compact object formation. In this paper, we present the first study of 
the persistent X-ray source Cygnus X-1 that takes into account of all available observational constraints. Our 
analysis accounts for three evolutionary phases: orbital evolution and motion through the Galactic potential 
after the formation of black hole (BH), and binary orbital dynamics at the time of core collapse. We find that 
the mass of the BH immediate progenitor is 15.0-20.0 Mq, and at the time of core collapse, the BH has 
potentially received a small kick velocity of < 77 km s~' at 95% confidence. If the BH progenitor mass is 
less than ^17 M©, a non zero natal kick velocity is required to explain the currently observed properties of 
Cygnus X-1. Since the BH has only accreted mass from its companion's stellar wind, the negligible amount 
of accreted mass is impossible to explain the observationally inferred BH spin of a, > 0.95, and the origin of 
this extreme BH spin must be connected to the BH formation itself. Right after the BH formation, we find 
that the BH companion is a 19.8-22.6 Mq main sequence star, orbiting the BH at a period of 4.7-5.2 days. 
Furthermore, recent observations show that the BH companion is currently super-synchronized. This super- 
synchronism indicates that the strength of tides exerted on the BH companion should be weaker by a factor of 
at least two compared to the usually adopted strength. 
Subject headings: binaries: close — X-rays: binaries — X-rays: individual (Cygnus X-1) 



1. INTRODUCTION 

In recent years, the number of observed black hole (BH) 
X-ray binaries (XRBs) has grown significantly. For these bi- 
naries, there exists a wealth of observation information about 
their current physical state: BH and donor masses, orbital pe- 
riod, donor's position on the H-R diagram and surface chem- 
ical composition, transient or persistent X-ray emission, and 
Roche lobe overflow (RLO) or wind-driven character of the 
mass transfer (MT) process. Furthermore, proper motions 
have been measured for a handful of these binaries (e.g. 
iMirabel et al.ll200Tl I200I [Mirabel & Rodriguesll2003h . To- 
gether with the earlier measurements of center-of-mass ra- 
dial velocities and distances, we can obtain information about 
the three-dimensional kinematic properties of these binaries. 
Given this plethora of observation results, the current ob- 
served sample of BH XRBs provides us with a unique oppor- 
tunity to understand the formation and evolution of BHs in bi- 
naries. This paper is the third in a series where we investigate 
in detail the BH formation in XRBs, especially focusing on 
the mass relationship between BHs and their immediate pro- 
genitors and the possible BH natal kick magnitude imparted 
during the core collapse event. 

In the first paper of this series, I Willems et al.l (I2005L here- 
after Paper I) showed how using the currently available con- 
straints one could uncover the evolution history of an XRB 
from the present state back to the time just prior to the core 
collapse event. They applied their analysis to the BH XRB 
GRO J 1 655-40. In the second paper. lFragos et all (120091 here- 
after Paper II) performed the same analysis for the case of the 
BH XRB XTE J1118H-480. In this work, we focus on the case 
of the BH XRB Cygnus X-1. The mass transfer mechanism 
in Cygnus X-1 is different from the XRBs studied in previous 



cases. Both donors in GRO J1655-40 and XTE J1118H-480 
are transferring mass to the BH under Roche lobe overflow, 
whereas the BH in Cygnus X-1 is accreting mass from the 
stellar wind of its companion. 

The plan of the paper is as follows. In Section 2, we review 
Cygnus X-l's currently available observational constraints. A 
general outline of the analysis used to reconstruct the system's 
evolutionary history is presented in Section 3, while individ- 
ual steps of the analysis are discussed in more detail in Section 
4-7. In Section 8, we derive constraints on the formation of 
the BH. The final section is devoted to a summary of our re- 
sults and discussion of some of the assumptions introduced in 
our analysis. 

2. observational constraints for cygnus x-1 

Cygnus X-1 was first detected in Aerobee surveys in 1964 
by iBowver et al.l (Il965h . Soon after the discovery, it was 
identified as an XRB, which cons isted of a compact object 
and a visible star HDE 226868 dMurdin & Webste^ [T97Tt 
I Webster & Murdinlll972l: lBoltonlll972ah . Spectroscopic ob- 
servations led lWalbornl(ll973h to classify H DE 226868 as an 
09. 7 lab supergiant. Bregman et al. ( 1973) estimated the dis- 
tance to be 2.5 kpc and set a lower limit of 1 kpc, based on the 
colors of field stars in the vicinity of the supergiant. Using a 
combination of data from David Dun lap Observatory (DDO) 
and the Royal greenwich Observatorv. lBoltoiil (1 1 972bl) derived 
the orbital period, eccentricity and systemic radial velocity 
(Vo) to be 5.5995 ± 0.0009 days, 0.09 ± 0.02 and -6.0 ± 0. 1 
km s"\ respectively. Based on the absence of X-ray and op- 
tical eclipses, the author gave a lower limit of 7.4 Mq on the 
mass of the compact object. This implied the compact object 
was too massive to be a white dwarf or a neutron star. Thus, 
the author confirmed iWebster & Murdin, (1972,) 's finding that 
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the compact object is a very strong BH candidate. 

Using the orbital period obtained from spectrometry and a 
range in th e assumed degree of Roche filling of the super- 
giant, |Gier& Bolton ( 1982, 1986a) found a lower mass limit 
of 7 Mq for the compact object. This confirmed that the com- 
pact object observed in Cygnus X-1 was a BH. The same 
authors also refined the orbital period and eccentricity to be 
5.59974 ±0.00008 and 0.021 ±0. 013, respective l y, and mea- 
sured Vq to be -2.0 ± 0.7 km s"'. iNinkov et all (119871) used 
the relationship between the equivalent width of the H7 spec- 
tral line and the absolute magnitude of early-type supergiants 
to estimate the d i stance as 2.5 ± 0.3 kpc . 

iHerrero et al.l (119951) performed a detailed spectroscopic 
analysis on the supergiant, and derived the masses to be 
10.1 M0 and 17.8 M© for the BH (Mbh) and the supergiant 
(M2), respectively, if an orbital inclination angle of 35° was 
assumed. Using the Isaac Newton telescope, LaSa la et aP 
(IT998.) measured the orbital period as 5.5997 ±0.0001 days 
and Vo as -5 .4 ± 0. 1 km s"' . With all the accumulated radial 
velocity measure ments and their own spectroscopy of the su- 
pergiant, [Broc^opp et al. ( 1999) refined the orbital period to 
5.599829 ±0.000'0'24. The proper motion of Cygnus X-1 was 
observed with the Very Large Baseli ne Interferometry ( VLBI) 
between 1988 and 2001 ( Lestrade eFal] IT999t iStirling et all 
120011: iMirabel & Rodriguesii2003i) . During this period, the 
system's position shifted at a rate of -4.2 ± 0.2 mas yr"' in 
right ascension (R.A.) and -7.6 ±0.2 mas yr"' in declination 
(dec). Meanwhile, a trigonometric parallax of 0.73 ± 0.30 
mas was also measured with VLBI, which gave a distance of 
1.4!°:^ kpc (ILestrade et al.l[T99 9l). 

By studying the spectra obtained with the 0.9 m coude 
feed telescope of Kitt Peak National Observatory, the 2.1 
m telescope of University of Texas McDonald Observatory, 
and the 1 .9m telescope of University of Tor onto David Dun- 
lap Observatory between 1998 and 2002, iGies et all (120031) 
derived Vq as -7 .0 ± 0.5 km s~^ and e s timate d Mbh/M2 
« 0.36 ± 0.05. Caballero-Nieve s et al.l (12009^ examined 
the supergiant' s ultraviolet spectra from the Hubble space 
telescope. Their results gave masses of 23;!;g and ll;!;^ 
M0 for the supergiant and the BH, respectively. On the 
other hand, Shaposhnikov & Titarchuk (2007) used the X-ray 
quasi-periodic oscillation and spectral index relationship and 
deduced Mbh to be 8.7 ±0.8 Mq, wh ich overlapped with the 
lower end of the Mbh range derived bv lCaballero-Nieves et al.l 
(l2009h . 

Recently, iReid et al] (1201 Ih measured the trigonometric 
parallax of Cygnus X-1 with the National Radio Astron- 
omy Observatory's Very Long Baseline Array (VLBA) and 
found a distance of 1.86;!;q [j kpc. The authors also reported 
proper motion measurements of Cygnus X-1, which were 
-3.78 ±0.06 mas yr"' in R.A and - 6.40±0.12 mas yr"' in 
dec. Meanwhile, Xiang etal] (1201 1) studied the X-ray dust 
scattering halo of Cygnus X-1 and determined the distance to 
be 1.81 ±0.09 kpc, after considering the compatibility with 
the parallax result. Bui lding on the trigonometric paral l ax dis- 
tance measurement of IReid etal.1 ( 120111) . lOrosz et aTl (120111) 
performed optical data modeling of Cygnus X- 1 , and found 
the mass of the supergiant to be 19.2 ± 1 .9 M© an d the black 
hole to be 14.8 ± 1.0 M©. Using the results o f iReid et al.l 
( 1201 Ih and lOrosz et all (120111) . iGou et al.1 ( 1201 Ih determined 
that Cygnus X- 1 hosts a near-extreme Kerr BH, with a spin 
parameter a, > 0.95. 

Unlike most of the XRBs known to host a BH, Cygnus X- 



1 is a persistent X-ray source. Sinc e the supergiant is cur - 
rently not overfilling its Roche lobe (iGies & Bolton|[T986al) . 
the observed X-rays are mainly powered by the accretion of 
stellar wind. The X-ray luminosity of Cygnus X-1 varies be- 
tween two discrete levels, namely the "hard (low) state" and 
the "soft (hi gh) state". As the system spends most of its time 
(-90%, see ICadolle Bel et al.l l2006) in the har d state, we fo- 
cus on the hard state X-ray luminosity (Lx). iFrontera et al.l 
(I2OOII) observed Cygnus X-1 with the Narrow Field Instru- 
ments of the BeppoSAX satellite at different epochs in 1996. 
The authors obtained the Lx (0.5-200keV) and the extrapo- 
lated bolometric luminosity (Lboi) as 2.0 x lO"'^ and 2.4 x 10-^^ 
erg s"', respectively, assuming a distance of 2 kpc. Using ob- 
servational data obtained by the Compton Gamma Ray Ob- 
servatory (CGRO) between 1991 and 2000. ' McConnell et all 
(1^02) derived Lboi to be ( 1 .62 - 1 .7( )) x lO" ergs"', with the 
distance to the source fixed at 2 kpc. ICadolle Bel et al.l (12006b 
observed Cygnus X-1 with the International Gamma-Ray As- 
trophysics Laboratory (INTEGRAL) between 2002 and 2004 
and measured Lx (20- 100 keV) as 6.5 x lO"*^ erg s"', as- 
suming a distance of 2.4 kpc. The authors also gave Lboi as 
2.2 X 10" ergs"'. 

For the systemic parameters relevant to our analysis, we 
adopt the most recent observational constraints, with the ex- 
ception of Lboi. We consider all the Lboi values mentioned 
above, assuming they represent the typical X-ray variabil- 
ity range for this system. After rescali ng their va lues to the 
parallax distance measurement by Reid et al.l (12011.) and con- 
sidering the uncertainty in that distance, we adopt Lboi to be 
(1.17-2.35)x 10^^ ergs"'. Foreaseof reference, our adopted 
observational constraints are summarized in Table 1 . 

3. OUTLINE OF ANALYSIS METHODOLOGY 

In our analysis, we assume that Cygnus X-1 formed in the 
Galactic disk, from the evolutio n of an isolated primordial bi- 
nary at solar metallicity. In fact. IMirabel & RodriguesI ( 120031) 
suggest that Cygnus X-1 belongs to Cygnus OB3 (Cyg OB3), 
which is an OB association located close to the Galactic plane. 
We also assume that no mass transfer via Roche lobe overflow 
occurred in the evolutionary history of this binary. 

According to our current understanding, in order to form 
a — 15 M0 stellar BH at solar metalicity, the BH progenitor 
in th e primordial binary nee ds to be more massive than 120 
M© (iBelczvnski et al.ll2010l) . Such a massive star loses its 
hydrogen rich envelope via stellar wind, and exposes its naked 
helium core. At the end of nuclear evolution, it collapses into 
a BH. During the core collapse event, the orbit is altered by the 
asymmetric mass loss from the system and a possible recoil 
kick imparted to the BH. If the binary survives through the 
core collapse event, angular momentum loss via gravitational 
radiation and tidal effects causes the orbit to shrink, although 
wind mass loss leads to orbital expansion. In the meantime, 
the more evolved BH companion is losing mass via its own 
stellar wind at a higher rate. The system becomes a BH XRB 
when the BH captures a non-negligible amount of mass from 
its companion's stellar wind. 

In this paper, we restrict ourselves to the formation of BH 
XRBs through the above evolutionary channel. Like the first 
two papers, our goal is to track the evolutionary history of 
Cygnus X-1 back to the time just prior to the core collapse 
event. Our analysis incorporates a number of calculations 
which can be summarized in four steps. 

First, we identify the current evolutionary stage of the BH 
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Table 1 

Properties of Cygnus X- 1 



Parameter 



Notation 



Value 



References 



Distance (kpc) 

Galactic longitude (deg) 

Galactic latitude (deg) 

Proper motion in R.A. (mas yr"' ) 

Proper motion in decl. (mas yr"') 

Systemic velocity (km s"') 

Orbital period (days) 

Orbital eccentricity 

Inclination angle 

Black hole mass (Mq) 

Black hole spin 

Companion mass (Mq) 

Companion Radius (Rq) 

Companion Luminosity (Lq ) 

Companion Effective temperature (K) 

Companion surface rotation speed (km s"') 

Bolometric luminosity of the X-ray source (erg s"') 



d 
1 
b 

IJ-RA. 

Vo 

Porb 

^orb 

(' 

Mbh 
a* 
M2 
R2 

u 

Tcff 
Vrot sin I 



1.86:«;if 

71.3 

+3.1 

-3.78 ±0.06 

-6.40 ±0.12 

-7.0 ±0.5 

5.599829 ±0.000016 

0.018 ±0.003 

27.06 ±0.76 

14.81 ±0.98 

> 0.95 

19.16±1.90 

16.50 ±0.84 

(1.91-2.75) X 10' 

30000 - 32000 

95 ±6 

(1-3-2-1)(t^)'x10" 



(9) 
(2) 
(2) 
(9) 
(9) 
(5) 
(1) 
(8) 
(8) 
(8) 
(10) 
(8) 
(8) 
(8) 
(8) 
(V) 
(3),(4),(6) 



References. — (1) Brocksopp et al. 1999, (2) Lestrade et al. 1999, (3) Frontera et al. 2001, (4) McConnell et al. 2002, (5) Gies et al. 2003, (6) Cadolle Bel et 
al. 2006, (7) Caballero-Nieves et al. 2009, (8) Orosz et al. 201 1, (9) Reid et al. 201 1, (10) Gou et al. 201 1 



companion, so that all the observational constraints are satis- 
fied. Under the assumption that the BH companion mass has 
not been altered by mass transfer in the past, we model it as 
an isolated star. Using a stellar evolution code, we calculate 
a grid of evolutionary sequences of isolated stars at differ- 
ent zero age main sequence (ZAMS) masses. We examine 
each sequence to find whether there exists a point in time that 
the calculated stellar properties, i.e. mass, radius, luminosity 
and effective temperature, are all simultaneously in agreement 
with the currently observed properties of the BH companion. 
If such a period of time exists, we classify that sequence as 
"successful". The current age of the BH companion can be es- 
timated from these successful sequences, and the time expired 
since the BH formation can then be derived by subtracting the 
approximate lifetime of the BH progenitor 

Next, we consider the kinematic evolutionary history of the 
XRB in the Galactic potential. Starti ng from the current loca- 
tion, we follow the methodology of iGualandris et al] ( l2005h 
and use the observed three-dimensional velocity to trace the 
Galactic motion of Cygnus X-1 backward in time. Together 
with the constraints on the current age of the system derived 
in the first step, this allows us to determine the location and 
velocity of the binary at the time of BH formation (we denote 
these as "birth" location and velocity). By subtracting the lo- 
cal Galactic rotational velocity at the "birth" location from 
the systems's center-of-mass velocity, we derive constraints 
on the peculiar velocity of the binary right after the formation 
oftheBH. 

In the third step, we analyze the orbital dynamics of the 
core collapse event due to mass loss and possible natal kicks 
imparted to the BH. In this paper, we refer to the instants 
right before and after the formation of the BH by the terms 
"pre-SN" and "post-SN", respectively. We start with the con- 
strained parameter space of (Mbh, M2) derived in the first step 
and perform a Monte Carlo simulation scanning over the pa- 
rameter space of the pre-SN binary properties. This parameter 
space is limited by requirements of orbital angular momentum 
and energy conservations, and by the post-SN binary peculiar 
velocity constraint derived in the second step. This calcu- 
lation yields a population of simulated post-SN binaries for 



each successful sequence. 

Finally, we follow the orbital evolution of these simulated 
binaries to the current epoch. Our calculation accounts for 
tides, wind mass loss, wind accretion onto the BH, and or- 
bital angular momentum loss via gravitational radiation. At 
the end of the calculations, we require agreement between the 
observed and calculated orbital period and eccentricity. 

4. MODELING THE BH COMPANION 

Under the assumption that the companion mass has not 
been altered by mass transfer in its past, we model the com- 
panion as an isolated star using a modified vers ion of the stel - 
lar evolution code EZ (originally developed b v lPaxtonll2004l) . 

We calculate the evolution of our stellar mo dels at solar 
metallicity, which is the same metallicity that ' Orosz et aTl 
(201 1) used in deriving the properties of the BH companion. 
When we place the companion's observational constraints on 
an H-R diagram, we find that the current location of the com- 
panion does not seem to be consistent with any evolutionary 
tracks calculated by the stellar evolution code. As shown in 
Figure 1, the companion is overluminous for a star of its mass. 
This cannot be explained by earlier mass transfer from the BH 
progenitor to the companion. Braun & Langer (1995) stud- 
ied the effects of mass accretion onto massive main sequence 
stars, and found that the accreting stars would not appear over- 
luminous for their new masses during the rest of their main 
sequence lifetime. If mass accretion leads to a so called "reju- 
venation" of the accreting star, which means its central hydro- 
gen abundance substantially increases, it would have the same 
luminosity as a star of its new mass. If rejuvenation does not 
occur, the accreting star would appear underluminous for its 
new mass during the rest of its main sequence lifetime. 

One possible solution for matching the observed compan- 
ion's luminosity is increasing the core overshooting parame- 
ter ttov to ^ 0.45. Although this value is relatively high, it 
is not unphysical. Claret (2007) compared the data from 13 
double-line eclipsing binary systems with theoretical predic- 
tions of stellar modeling and found aov could be as high as 
0.6 for massive stars. We vary Oov from 0.35 to 0.5, in steps 
of 0.01. We note that the need for such higher values of aov 
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Figure 1. The evolutionary tracks for isolated stars on the H-R diagram. On 
each track, the mass of the star in M© is indicated at various points. The 
gray shaped area represents the observational constraints of the BH compan- 
ion. At Teff ~ 3100()K. the model with an initial mass of 20 Mq (dashed 
line) would have a mass Rd 19.64 Mq, which is in good agreement with the 
measured mass of the companion. However, it does not match the measured 
luminosity. On the other hand, the model with an initial mass of 30 Mq (dot- 
ted line) match the measured luminosity at Tjff Si 31000 K, but the mass of 
the star does not match the companion's. The model with an initial mass of 
20 Mq and cJov = 0.45 (solid line) could match both the measured mass and 
luminosity of the companion at T^ff Ri 31000 K. 

in the modeling of massive stars may very well be connected 
to the significant presence of internal rotation and associated 
rotational mixing. Effectively increasing Oov leads to stronger 
internal mixing and in a way allows the stellar model to be- 
have more like a rotating model. 

Besides the observational constraints on the companion's 
properties, there are three additional constraints. The first one 
comes from the fact th at the companion is cur rently not over- 
filling its Roche lobe (iGies & Boltonlll986ah . Thus, we re- 
quire the stellar radius R2 in our models to be 



R2 < Aoi-brEgg + A/; 



(1) 



where Tess is the effective Roche lobe radius given by 
lEggletod (11983). Here, we make an approximation that the 
orbit is circular and synchronized. The parameter AR is a 
constant accounting for the difference in the calculated stel - 
lar radii among stellar evolution codes (IValsecchi et al.l2()10l) . 
Weset ARto2.5R0. 

Another constraint is that the calculated bolometric lumi- 
nosity (Lboi) resulting from the stellar wind accretion process 
needs to fall within the observational range, which is ( 1.17- 
2.35) X 10" erg s"'. By adopti ng thelBondi & Hovld (11944 
accretion model and following iBelczynski et al. ( 20081) , the 
orbital-averaged accretion rate is given by 






\ wind 



^wind 



M2 



2AL(l+^')'/' 



(2) 



Here, Mbh is the BH mass in our models, which varies within 
the Icr range of the observational constraint, in steps of 0.098 
M0. Since the total mass that the BH could have accreted 
from its companion stellar wind is negligible, Mbh in each 
evolutionary sequence is fixed throughout our analysis. M2 is 



the wind mass loss rate of the companion in our models. Fwind 
is a parameter such that M^cc never exceeds -O.8M2, and awind 
is the accretion efficiency, which varies between 1.5 and 2.0 
(JBoffin & Jorissen 1988). Aorb and Coib are the orbital semi- 
major axis and eccentricity, respectively. Aorb is derived from 
the mean measured orbital period Porb, which is 



iorh ■ 



G(MBH+M2)Plh 



47r2 



(3) 
Corb is set 

* wind ^^" 



where M2 is the companion mass in our models. 

equal to the mean measured orbital eccentricity. 

notes the wind velocity. V^ equals to V^jj/V^ji^^, where V|jj 

is the orbital velocity square of the BH and is approximated 

as G(Mbh + M2)/Aoi-b. We adopt the spherically symmetric 

wind velocity law given iu iLamers & Cassinelli ( 1999), 



V„i„d(r) = Vesc + (Vex 



Vesc) ( 1 - y 



/3 



(4) 



where r is the distance from the companion to the BH and 
is set equal to Aorb- /3 is a free parameter varying from 0.6 
to 1.6 (Gies & Bolton 1986b; Lamers & Leitherer 1993), in 
steps of 0.1. Voo is the wind velocity at infinity, while Vesc 
is the effective escape velocity at the surface of the compan- 
ion. Within the typical r ange of O star surface t emperature, 
Voo is s caled as 2.65Vesc (i Kudritzki & Pulsll2000l) . Following 
iLamers & Cassin elli ( 199^ 



where 



Vesc = ^/2(l-^e)GM2/R2: 
p ^ CreL2 

AttcGM2 



(5) 
(6) 



is the mass correcting factor for the radiative force due to 
electron scattering, and c is the speed of light in vacuum. 
ILamers & Leitheren (119931) scaled the electron scattering co- 
efficient per unit mass ag as 



cr, = 0.401 



1 + qe 
1 + 3e 



(7) 



where q is the fraction of He"^ and (1 - q) is the fraction of 
He+, with q = 1 if Tett > 35,000 K, q = 0.5 if 30,000 K < Teff 
< 35,000 K, and q = if Teff < 30,000 K. The abundance ratio 
e = He/(H -1- He) is fixed at 0.15, which is appropriate for an 
O star with a spect ral type of Class I . Using Mace from equa- 
tion (2), we follow Be Iczynski et al .' (2008) and calculate the 
bolometric luminosity resulting from the companion's stellar 
wind being accreted onto the BH as 



i'bol = 



1 GMBHMa, 



(8) 



where R^cc denotes the radius of the accretor. For the case 

of BH, Race is the radius of the inner most stable circular or- 

bit, wh ich we calculate with Equation (2.21) in lBardeen et aTl 

19721). Given t he observationally inferred spin a* > 0.95 

Gouetal.ll2011h . we adopt the median a* = 0.97 and find 



-'^acc — 2.57 



Mr 



Mp 







km. 



(9) 



This calculated Lboi needs to fall within the observational 
range. 
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Figure!. Systemic behavior of two selected evolutionaiy sequences, wliicli have the same cjov = 0.44, Mbh = 14.81 Mq, a„i„^ = 1.5, and /3 = 1.0. Sequence 
1 (solid) and 2 (dashed) have M2,zanis of 21 and 22 Mq, respectively. The left panel shows the evolutionary tracks on the H-R diagram, while the middle panel 
illustrates the behaviors of the mass and the radius of the star. The right panel shows the variations of the calculated stellar luminosity and /l, where /l is defined 
in Equation (10). The gray shaded areas represent the observational constraints on the relevant quantities, and the thick part of the evolutionary tracks indicates 
the part of the sequence that the observational constraints on the H-R diagram are satisfied. 

The last additional constraint is that the observational con- 
straints on L2 and L[,oi have to be evaluated at the same distant 
estimation. To examine this, we calculate the ratio 



0.46 



fi = 



/ L2 \ / Lboi 



\10^LqJ \10^''ergs- 



(10) 



which is independent of distance. From Figure 1 in 
IOroszetal.l(l20nl) . L2 is 2.09 x 10^ Lq at Teff = 30000K, and 
is 2.51 X 10^ L0 at Teff = 32000K, assuming a distance of 1.86 
kpc. Together with the measured range of Lboi rescaled at the 
same distance estimation, the upper and lower limits of /l are 
1.01 and 1.90, respectively. We can assure that both luminos- 
ity constraints are evaluated at the same distance estimation if 
/l falls within that range. 

In order to find the current evolutionary stage of the BH 
companion, we apply these constraints to a set of evolution- 
ary sequences, which cover the parameter space of the com- 
panion's ZAMS mass (M2,zams), <^ov, Mbh, awind, and /3. For 
each sequence, we find whether there exists a point in time 
that the calculated properties simultaneously satisfy all ob- 
servational constraints: the BH companion's mass, luminos- 
ity, temperature, and radius, Lboi, /l, and not overfilling the 
Roche lobe of the BH companion. Similar to the Roche lobe 
constraint, we also consider an uncertainty of ±2.5 R© in the 
calculated stellar radii when we apply the observational con- 
straint of the BH companion's radius. If such a period of time 
exists, we classify that evolutionary sequence as "successful". 
The behavior of some relevant parameters is illustrated in Fig- 
ure 2 for two selected successful sequences that are chosen 
mainly to provide a clear and instructive picture. The dis- 
played sequences therefore do not represent our best possible 
matches to the observed properties of Cygnus X-1. Figure 3 
shows the parameter space of Ma.zams, ctov, and Mbh covered 
by all successful sequences. For Owind and /3, the successful 
sequences covered the entire allowed parameter space, which 
are 1.5 < Owind < 2.0 and 0.6 < ^ < 1.6. 

The current age of the BH companion could be derived 
from the time interval at which all observational constraints 
are satisfied. Assuming that the BH progenitor and its com- 
panion formed at the same time, we could compute the time 
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since the BH formation (tjys) by 



20 21 22 23 

M2,zams (^o) 
IV, and Mbh covered by all 



tsys - h - fBH , 



(11) 



where tBH i s the approxima t e lifet ime of the BH progenitor 
We follow Belczvnski et al.l (120101) to calculate Mbh and tBH 
for different progenitors using the stellar evolution code SSE 
(Hurley et al. 2000), and adopting the mass loss prescriptions 
which were classified as "Vink et al. Winds". The calculated 
tBH are fit as a function of Mbh, 



ffiH 



Me 



■ + 3.341, 



lOVs 19.26-4.902 f^") +0.3841 (^ 

(12) 
for Mbh > 9.5 Mq. Figure 4 shows the variations of X.2 and 
tBH against M2.zams- We find that tsys is between 4.8 and 7.6 
Myr. 
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Figure 4. The variations of t2 (circles) against Mizams- The gray shaded 
region indicates the range of ten for the corresponding successful sequences, 
which is calculated by Equation (12). The difference between t2 and tBH 
gives tsys. 



5. KINEMATIC HISTORY IN THE GALAXY 

Here, we assume that Cygnus X-1 formed in the Galactic 
disk. The consideration of Cyg OB3 being the parent associa- 
tion of Cygnus X-1 is discussed in Section 9.2. Given the ob- 
served position and measured proper motion of Cygnus X-1, 
we derive the post-SN peculiar velocity of the binary's center- 
of-mass by tracing its orbit in the Galaxy back to the time of 
BH formation. We describe the motion of the binary with re- 
spect to a right-hand Cartesian reference frame, whose origin 
coincides with the Galactic center. The Z axis points to the 
northern Galactic pole, while the X axis points in the direc- 
tion from the projected position of the Sun onto the Galactic 
plane to the Galactic center. In this reference frame, the Sun 
is located at (Xp. Yp,, Z^;,) = (-8.5,0 0.03) kpc jJ oshi 2007; 
iGhez et aLllIOOl lGillessen"erani200^ iReid et alJl2009l) . and 
has a peculiar motion (Uq, V q, W©) = (11.1, 12.24, 7.25) 
km s~' dSchonrich et al.ll20l"oh . Cygnus X-1 is currently lo- 
cated at a distance of 1.86;^}j kpc from the Sun, with a 
Galac tic longitude / = 71 . 3°, and a Galact ic latitude b = 
3.1° d Lestrade et all [19991: IReid et all lIoTli) . This means 
Cygnus X-1 is currently '^ 130 kpc above the Galactic plane. 

To model the Galaxy, we adopt the Galactic potential of 



Carlberg & Innaneri (11987) with updated model parameters of 



Kuijken & Gilmore ( 1989). The equations governing the sys- 



tem's motion in the Galaxy are integrated backward in time, 
up to the time corresponding to the current system's age tsys 
given by the successful sequences. We follow the method- 
ology of G ualandris et al.. (2005) to initialize the parameters 
for the integration, which accounts for the uncertainties in the 
estimated distance and measured velocity components. We 
generate the initial system's position by the Galactic coordi- 
nates (/, b) and a random distance drawing from a Gaussian 
distribution. We generate initial system's velocity by drawing 
randomly the proper motions (/iR.A.,Mdeci.) and heliocentric 
radial velocity (Vq) from Gaussian distributions. The current 
system's age is uniformly distributed between 4.8 and 7.6 Myr 
(see Figure 4). 
Figure 5 shows the possible positions of Cygnus X-1 at the 
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Figure 5. Upper panels: The grey dots illustrate the possible locations of 
Cygnus X-1 at the birth time of the BH, obtained from 3,000 integrations of 
its trajectory backwards in time. The initial conditions of the integrations are 
generated randomly using the methodology described in Section 5. The plus 
signs indicate the cuirent location of Cygnus X-1, derived from the mean 
distance of 1.86 kpc. The crosses represent the current location of Cyg OB3 
center, with an adopted distance of 2 kpc. Lower panel: The distribution of 
post-SN peculiar velocities Vp^;. pojtSN against the time expired since the BH 
formation. 

time of BH formation, obtained from integrating 3,000 trajec- 
tories backwards in time. As there is no trajectory crossing the 
Galactic plane and the end points of all trajectories fall within 
1 10 pc from the Galactic plane, we consider each end point as 
a possible birth site of the BH. The post-SN peculiar velocity 
VpecpostSN of the binary is obtained by subtracting the local 
Galactic rotational velocity from the center-of-mass velocity 
of the binary at the birth sites. We find Vpec.postSN ranges from 
22 to 32 km s"' and is time independent. The distribution of 
VpecpcstSN against the time expired since the formation of the 
BH are displayed in Figure 5. 

6. ORBITAL DYNAMICS AT CORE COLLAPSE 

For each of the successful sequence, we perform a Monte 
Carlo simulation which consists of twenty million pre-SN bi- 
naries. The properties of the BH progenitor's companion are 
taken from the stellar model of that sequence, at the time when 
the age of the star is equal to Ibh- During a supernova (SN) ex- 
plosion, the mass loss from the system and possibly the kick 
imparted to the BH change the binary's orbital parameters. 
The pre- and post-SN component masses, orbital semi-major 
axis, and orbital eccentricity are related by the conservation 
laws of the orbital energy and angular momentum. In the fol- 
lowings, we add the subscripts "preSN" and "postSN" to the 
notations of the orbital elements to distinguish between their 
values just prior and right after the SN explosion that formed 
the BH. 

We start with seven free parameters: the BH immediate 
(He-rich) progenitor mass (Mhb), pre-SN orbital semi-major 
axis (ApreSN) and eccentricity (CpreSN), the mean anomaly (m), 
the magnitude (V^) and direction {9, (p) of the kick veloc- 
ity imparted to the BH. 9 is the polar angle of the kick with 
respect to the relative orbital velocity of the BH progenitor 
just prior to the SN explosion, and (f> is the corresponding az- 
imuthal angle (see Figure 1 in .Kalog era i 200 0., for a graphical 
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representation). The first five parameters are drawn from uni- 
form distributions, while the last two are drawn from isotropic 
distributions. It is obvious that the progenitor must of course 
be more massive than the BH, but there is no absolute upper 
limit for the progenitor mass. We adopt Mhe <20 Mq, and 
provide a discussion on this upper limit in Section 9.1. 

The relations between pr e- and post-SN parameters have 
been derived bv lHillsl (Il983h : 



n + ^He.preSN + 2yiyHe,preSN COS 6* 

= G(Mbh+M2) 



, '" ApostSN / 
G(MBH+M2)ApostSN(l -SpostSN) 

= r^\yls\n^ecOS^(t>+ [sinf/'(VHe,preSN + V4COS 

- Vk COS il) sin 9 sin (j) 



(13) 



(14) 



Here, r is the orbital separation between the BH progenitor 
and its companion at the time of SN explosion, 

'" = Ap,eSN(l - SpreSNCOS -EpieSN), (15) 

where E is the eccentric anomaly, and is related to m as 

m = E - e smE. (16) 

Vne.preSN IS the relative pre-SN orbital velocity of the BH pro- 
genitor. 



Vi 



He.preSN ■ 



G(MHe+M2)(---^ 
r ApreSN 



ll/2 



(17) 



The angle ip is the polar angle of the position vector of the BH 
with respect to its pre-SN orbital velocity in the companion's 
frame. It is related to the pre-SN parameters as 



'•^^He^preSNSinV 



G(MHe+M2)Ap,eSN(l-4eSN)- (18) 

Since the core collapse is instantaneous, r remains unchanged. 
This gives a constraint 



r = ApreSN(l - epreSNCOSiipieSN) 
= ApostSN(l - epostSNCOsZipostSN), 



(19) 



which needs to be satisfied with | cosZSpostSN I =$ 1 • 

The mass loss from the system and a natal kick imparted 
to the BH can induce a post-SN peculiar velocity (Vpec,postSN) 
at the binary's center of mass. Its magnitude is determined 
by following Equations (28)-(32) in Paper I, and is required 
to fall within the range derived in Section 5, which is 22-32 
km/s. 

In addition, there are two more restrictions on the proper- 
ties of pre- and post-SN binary components. First, we require 
that both components have to fit within their pre- and post-SN 
Roche lobe at periapsis. We impose this condition to avoid 
complications arising from mass transfer induced changes in 
the stellar structure of the MS companion, that later becomes 
the BH companion of the XRB. To calculate the Roche lobe 
radius of each component in eccentric pre- and post-SN or- 
bits, we adopt the fitting formulae of Sepinsky et al. (2007). 
When calculating the pre-SN Roche lobe radii, we assume 
that the pre-SN orbit is pseudo-synchronized. Again, due to 
the difference in calculated stellar radii among stellar evolu- 
tion codes, we consider an uncertainty of ±2.5 R© on the 
companion radius (.Valsecchi et a l. 2010,). The radius of the 



BH i mmediate progenito r can b e approximated by Equations 
(3) in iFrver & Kalogeral ( 119971) . since we assume that it is 
a Helium star Second, the pre-SN spin of the BH imme- 
diate progenitor and its companion need to be less than the 
breakup angular velocity fie « (GM/R^)'/-. As the calculated 
stellar radius R? assoc iates with an uncertainty AR = 2.5R0 
dValsecchi et al.ll20Tol) . 



IGM7 



3AR 



(20) 



^r 



for the BH companion. 



7. ORBITAL EVOLUTION AFTER THE SN EXPLOSION 

The orbital evolution of the simulated binaries, which are 
generated from the Monte Carlo simulations described in Sec- 
tion 6, is calculated up to the current epoch. After the forma- 
tion of the BH, the orbital parameters of the binary are subject 
to secular changes due to the tidal torque exerted by the BH 
on its companion, and due to the loss of orbital angular mo- 
mentum via gravitational radiation and stellar wind. Since 
the tidal interactions depend on both the orbital and rotational 
properties of the MS companion, the star's rotational angular 
velocity (fi) right after SN explosion that formed the BH en- 
ters the problem as an additional unknown quantity. Here we 
assume the rotational angular velocity of the BH companion is 
unaffected by the SN explosion, and is pseudo-synchronized 
to the pre-SN orbital frequency. The system of equations gov- 
erning the tidal evolution of the orbital semi-major axis A, 
eccentricity e, and the BH co mpanion's rotational angular ve- 
locity VL has been derived by lHud (Il981h : 



dA 
dt 



tides 



_^MbhMbh + M2 [Ri 

'r M2 M2 

A 



(l-e2)i5/2 



1\ /1 2\3/2 ^ /- 2^ ^^ 



f,[e^)-{l-ey''h{e^)- 



(21) 



de\ ^ _ 27 ^2 Mbh Mbh +M2 f R2 



dt 



dn 

'df 



tides 



T M2 M2 



(1-^2)13/2 



A 
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hi^')-^A^-^'f'f4^') 
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n 
(22) 



= 3- 



kl fMBH\ M2RI (R2 



tides 



M2 



h 



(l-e2)' 



2\b 



2^3/2, 



f,{e')-{l-eyy,{e') 



n 



(23) 



Here, k2 and I2 are the apsidal-motion constant and moment 
of inertia of the MS companion, respectively. T is a charac- 
teristics timescale for the orbital evolution due to tides, and n 
= 27r/Porb is the mean orbital angular velocity. The coefficient 
fu nctions /} ( e^) for i = i, 2, . . . , 5 are given in Equations (11) 
in lHutl ( ir981[) . As the BH companion in Cygnus X-1 is a mas- 
sive MS star with a radiative envelope, the factor k2/T can be 
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approximated as 



= 1.9782 X lO"* 



rad 






Rq 

A 



5/2 



M2 V^^ fMBH+M2^^^^ 



El yr 



(24) 



^MqJ \ M2 
The constant E2 comes from a fit to the tables in lClareti (120041) 



logEi 



t/t„ 



2.20489- 1.89579(f/f„„) 



-5.51039, 



(25) 



for 15.85 < M2,zams < 25.12 M0. Here, t^s is the main se- 
quence lifetime. We define the end of the main sequence as 
the hydrogen abundance at the core being less than 0.01. 

To follow the secular changes of the orbital parameters 
associated with emissions of gravitational waves, w e adopt 
Equations (35) and (36) in 'J unker & Schaferi (119921) . which 
are derived up to 3.5 post-Newtonian order 

The rates of change in A and e due to wind mass loss 
and wind accretion onto th e BH are determine d by following 
Equations (15) and (16) in (iHurlev et al.ll2002l) . 



— =-A 
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dt, 


Wind 


'de^ 


wind 



M2 



1+e^ 



Mbh + Mi V^bh MBH + M2/ 1-e 



M„ 



2 
(26) 



-eMa, 



1 



1 



MBH+M2 2Mbh 



(27) 



The mass loss via stellar wind also induces a loss in the 
spin angular momentum of the BH companion. iHurlev et al.l 
(l200ft) showed that if all the mass is lost uniformly from a 
thin shell at the surface of the MS star. 



^pm =Jj^h^) = \M2Rln, 



(28) 



where 72,spin is the spin angular momentum of the BH com- 
panion. 

For each of simulated binaries, we follow the secular 
changes of its orbital properties due to all the mechanisms 
mentioned in this section. The properties of binary com- 
ponents are adopted from the corresponding successful se- 
quence. Unlike finding Vpec.postSN in Section 5, the orbital 
evolution of the binary goes forward in time, from ^bh to ti- 
Within this period of time, the BH companion has to always 
fit within its Roche lobe at periapsis. In other words, its cal- 
culated radius is constrained to be less than the R oche lobe 
radius at periapsis given by iSepinskv et alj ( 120071) . Again, 
we allow an uncertainty of ±2.5 R©, due to the differ- 
ence in calculated ste llar radii among stellar evolution codes 
dValsecchi et al.ll20r0 ). Furthermore, the rotational angular 
velocity of the BH companion has to be smaller than the 
breakup angular velocity Q,c. If the orbital period and eccen- 
tricity of the simulated binary at f2 match the measured values 
of Cygnus X-1, we classify that binary as a "winning binary". 
Figure 6 illustrates the time evolution of orbital parameters 
for one selected winning binary, and shows that the change in 
the semi-major axis is mainly determined by the stellar wind 
mass loss, while the change in the eccentricity is overwhelm- 
ingly dominated by the tidal effects. 

8. PROGENITOR CONSTRAINTS 
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Figure 6. The orbital evolution of a selected winning binary. Right after 
the formation of the BH (t = 3.8 Myr), this binary consists of a 14.8 M© BH 
and a 2 1 .7 M0 main sequence star. The top panels show the time evolution of 
orbital period and eccentricity. The bottom panels show the rate of changes of 
the semi-major axis A and eccentricity e due to tidal effects (solid line), wind 
mass loss and wind accretion onto the BH (dashed line), and gravitational 
radiation (dotted line). 



The elements presented in the previous sessions can now be 
combined to establish a complete picture of the evolution of 
Cygnus X-1 and the dynamics involved in the core collapse 
event that formed the BH. After finding the successful evo- 
lutionary sequences that satisfy all the observed properties of 
the BH companion and the bolometric X-ray luminosity as 
discussed in Section 4, we trace the motion of the system in 
the Galaxy back in time to the forma t ion of the BH. We adopt 
the methodology of lGualandris et all ( l2005h to account for the 
uncertainties in the measured distance and velocity compo- 
nents of Cygnus X-1. The time of BH formation is different 
for each successful sequence. It is estimated by the BH mass 
of the sequence, which connects to an approximate lifetime 
of the corresponding BH progenitor This procedure gives us 
a constraint on the system's peculiar velocity right after the 
BH formation. We then perform Monte Carlo simulations on 
the orbital dynamics at core collapse for each successful se- 
quence. There are seven free parameters: the BH immediate 
progenitor mass, the pre-SN orbital semi-major axis and ec- 
centricity, the mean anomaly, the magnitude of kick veloc- 
ity imparted to the BH, and the two angles specifying the 
direction of the kick velocity. The Monte Carlo simulations 
produce a population of simulated binaries, which satisfy the 
post-SN system's peculiar velocity constraint derived already. 
Last, we evolve the orbits of these simulated binaries forward 
in time to the current epoch. If the orbital period and eccen- 
tricity of the simulated binary at current epoch match the mea- 
sured values of Cygnus X-1, we classify that simulated binary 
as a "winning binary". The results presented in what follows 
are derived from the winning binaries of all successful se- 
quences. 

In Figure 7, we present the probability distribution func- 
tions (PDFs) of the BH immediate (He-rich) progenitor mass 
(Mhc) and natal kick magnitude (Vk). We find Mhc to be in 
a range of 15.0-20.0 M©, and Vk to be < 77 km s"^ both 
at 95% confidence. Figure 8 illustrates the 2D joint Vk-MHe 
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Figure 7. The probability distribution functions of tlie BH immediate (He- 
rich) progenitor mass (Mhc) and natal kick magnitude (V|[) imparted to the 
BH. 
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Figure 8. The 2D joint Vk-MHc confidence levels: 68.3% (red), 95.4% (yel- 
low), and 99.7% (blue). 

confidence levels, which shows that if Mhc is less than ^^17 
M0, the BH might have received a non-zero natal kick at the 
core collapse event. For small Mhe, a minimum V^ of ^ 55 
km s~' is necessary for explaining the current observed prop- 
erties of Cygnus X- 1 . Furthermore, both the Mhc PDF and the 
2D joint Vk-Mne confidence levels show that the maximum 
Mhc is constrained by our adopted upper limit of 20 Mq. We 
impose this limit based on the physics involved in the evolu- 
tion of massive stars. A discussion on this limit can be found 
in Section 9.1. Given our understanding of mass loss from 
Helium stars, it seems that the BH has potentially received a 
small natal kick velocity of < 77 km s"' (95% confidence) 
during the core collapse event. 

Based on the dynamical model of lOrosz et aTl (1201 Ih . 
iGou et all (1201 li) found that the BH in Cygnus X-1 has a spin 
parameter a* > 0.95 at Sc To determine whether the BH was 
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Figure 9. The 2D joint ApieSN-SpreSN confidence levels: 68.3% (red), 95.4% 
(yellow), and 99.7% (blue). 

born with an extreme spin, we first estimate how much mass 
the BH could have accreted from its companion's stellar wind 
since the time of BH formation. The winning binaries of all 
successful sequences show that at maximum the BH has ac- 
creted ~ 2 X 10"^ M0. Since it is impossible to spin the BH up 
to a^, > 0.95 by accreting that negligible amount of mass, the 
BH needs to have an extreme spin at birth. This high spin has 
implications a bout BH forma t ion an d the role of rotation in 
core collapse. lAxelsson et al.l (120111) also concluded that the 
observed spin connects to processes involved in core collapse, 
and is not likely to originate from the synchronous rotation of 
the BH progenitor. 

Besides the constraints on the BH formation, our results 
also shed light on the evolutionary picture of Cygnus X- 1 . We 
find that right after the formation of the BH, the BH compan- 
ion has a mass of 19.8-22.6 M©, in an orbit with period of 
4.7-5.2 days. Since then, the orbital separation of Cygnus X- 
1 has been increasing with time, as the rate of change in the 
semi-major axis is dominated by the influence of stellar wind 
mass loss from the system. On the other hand, the orbital ec- 
centricity has decreased slightly since the BH formation. This 
is because the tides exerted on the companion by the BH, as 
the dominant mechanism of circularizing the orbit, are not 
strong enough to decrease the orbital eccentricity significantly 
within the time period of several million years since the time 
of BH formation. We find that CpostSN ranges from 0.015 to 
0.022. However, this does not suggest that CpreSN has to be 
small. An eccentric pre-SN orbital could become fairly cir- 
cular if there is a natal kick imparted to the BH at the right 
direction. As illustrated in Figure 9, there are winning bina- 
ries with CpreSN being as high as ^ 0.53. 

9. SUMMARY & DISCUSSION 

In this paper we constrained the progenitor properties and 
the formation of the BH in the persistent XRB Cygnus X- 
1. Our analysis accounts for the orbital evolution and mo- 
tion through the Galactic potential right after the BH forma- 
tion, and the binary orbital dynamics at the time of core col- 
lapse. We find that the mass of the BH immediate progenitor 
falls within a range of 15.0-20.0 M0 at 95% confidence. 
We note that the maximum progenitor mass is constrained by 
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our adopted upper limit, which is discussed in Section 9.1. 
The BH has potentially received a small natal kick velocity 
of < 77 km s"' at 95% confidence. In fact if the progenitor 
mass is less than '-^ 17 M0, a non zero natal kick velocity 
is necessary to explain the currently observed properties of 
Cygnus X-1. Since the BH has only accreted mass from its 
companion's stellar wind, the total amount of mass accreted 
since the BH formation is less than ■^ 2 x 10"^ M0. This in- 
dicates that the o bservationally inferred BH spin of a, > 0.95 
JGouetal .1120 111) cannot be explained by mass accretion and 
has to be natal. This high spin has implications about BH for- 
mation and the role of rotation in core collapse. Right after the 
BH formation, the BH companion has a mass of 19.8-22.6 
M0, in an orbit with period of 4.7-5.2 days and eccentricity 
of 0.015-0.022. Although the post-SN orbital eccentricity 
is small, the pre-SN orbit can potentially be fairly eccentric. 
This is possible if the BH receives a natal kick velocity at the 
right magnitude and direction. 

The formation of the BH in Cygnus X- 1 has been 
previously studied by | Nelemans et aP (Il999h and 
iMirabel & RodriguesI (|2003|K Both studies assumed 

symmetric mass loss during the core collapse event, and 
considered only the binary orbital dynamics at the time of 
core collapse. Comparing with these two earlier studies, 
we consider the possible asymmetries developed during the 
core collapse event and the evolution of the binary since the 
BH formation. It is important to note that these two earlier 
studies do not consider the multitude of the observational 
constraints taken into account here and hence the suggested 
progenitors are not complete solutions for the evolutionary 
history of Cygnus X- 1 . 

Finally, we discuss some of the assumptions introduced in 
our analysis in the following sub-sections. 

9.1. Maximum BH Progenitor Mass 

Unlike the case of GRO J1655-40 studied in Paper I, the 
analysis of orbital dynamics during the core collapse event 
does not give an upper limit on Mm- Instead, we have con- 
servatively adopted an upper limit of Mee <20 Mq, based on 
physics involved in the evolution of massive stars. As men- 
tioned in Section 6.1 of Paper 11, by evolving a ZAMS star 
of '--^ 100 Mq at solar metallicity, the maximum Helium star 
mass one can achieve is ^15 M0 when including moderate 
stellar rotation, and ^17.5 Mq when assuming no stellar ro- 
tation. When adopting the upper limit of 17.5 Mq, the lower 
limit of Mhb decreases slightly to 14.6 Mq and the range of 
Vk becomes 14-81 km s~', both with 95% confidence. This 
range of V^ still suggests that the BH in Cygnus X-1 received 
a low kick during the core collapse event. 

9.2. Association with Cyg 0B3 

The center of Cyg OB3 locates at I = 72.8° and b = 
2.0°, and at a distance of 1.4- 2.7 kpc away from the 
SunjMassev et al. 1995; Dambis e t al. ll2001l;lMel'Nik et aTl 
l200iriMerNik & Dambis ii2009i) . When comparing that to 
the location of Cygnus X-1 (Table 1), it is clear that not 
only their Galactic coordinates are close to each other, but 
also their distance estimations overlap with each other Fur- 
thermore, the measurements of proper motion and radial ve- 
locity show that Cygnus X-1 is moving as the members of 
Cyg OB3 (Dambis et al. ^OOlt IMirabel & Rodrigue"sll200l 
Mel'Nik & Dambis 2009V Based on these observations, 
Mirabel & Rodrigu es C2003,) argue that Cyg OB3 is the par- 
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Figure 10. The same plot of 2D joint Vi^-MHe confidence levels as Figure 8, 
calculated with Vpec.postSN < 10 km s"' instead of the original range derived 
Section 5. 

ent association of Cygnus X-1. This infers that Vpec,postSN due 
to the core collapse event has to be small. If we change the 
constraint on Vpec.postSN to < 10 km s"', we find Mhc to be in 
a range of 13.9- 16.9 Mq and Vk to be < 24 km/s, both at 
95% confidence. Besides the change in 95% limits, non-zero 
BH natal kicks are not needed for progenitors of Mhc < 17 
Mq in order to explain the observed properties of Cygnus X- 
1, but become necessary for Mee > 17.5 Mq (see Figure 10). 
Also, we note that a relatively small change on the range of 
VpecpostSN affects the derived constraint on Vk qualitatively. 

9.3. Super-Synchronized Orbit 

After considering several previous measurements of 
the BH companion's surfa ce rotation speed (Vrotsin/), 
'Caballero-Nieves et all (120091) adopted Vi-ot sin / = 95 ± 6 km 
s~'. Orosz et al. (20111) found that the ratio of the BH com- 
panion's spinning frequency to the orbital frequency (/j^) was 
1.400 ±0.084, which was derived based on their results of the 
inclination angle / = 27°. 06 ±0°. 76 and the companion ra- 
dius R2 = 16.5 ± 0.84. This indicates that the BH companion 
is super-synchronized. We note that with the analysis pre- 
sented here, we find none of our winning binaries have super- 
synchronized BH companions at the current epoch. They are 
all sub-synchronized with fn reaching ^ 0.87 at maximum. 

In an effort to examine how our standard assumptions can 
be modified and investigate whether super-synchronism is at 
all allowed by the models as indicated by the observations, 
we make two modifications to our analysis. We first remove 
the assumption that the pre-SN orbit is pseudo-synchronized, 
and randomly distribute the pre-SN spin of the BH companion 
between zero and its breakup angular frequency fl^. Next, we 
reduce the secular changes of the orbital parameters due to 
the influence of tides by multiplying the right hand side of 
Equation (21) - (23) by a constant /tide. 

As shown in Figure 11, by allowing the pre-SN spin of 
companion to be greater than pseudo-synchronization and 
keeping the tidal strength unchanged (i.e. /tide = 1 0), the max- 
imum fn of the winning binaries increases to ^ 1.2. Although 
it is getting close, this value is still below the observation- 
ally inferred one. Together with a weakened tidal strength of 
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Figure 11. The probability distribution functions of the ratio (/q) of 
the companion's spinning frequency to the orbital frequency at cun'ent 
epoch, which are calculated with models of fm^ = 0.2 (solid/blue), 0.5 
(dashed/green), and 1 .0 (dotted/red). All three set of models have the pre-SN 
spin of the companion randoinly distributed between zero and the breakup 
angular frequency f2c . 

/tide = 0.2 and 0.5, we comfortably find winning binaries with 
/si = 1 .4. Furthermore, the minimum pre-SN surface rotation 
speed of the companion in those winning binaries are ~ 500 
and 700 km s"' for /ide = 0.2 and 0.5, respectively. Given 
the uncertainties in the physics of fast rotating massive stars, 
it seems that super-synchronism in Cygnus X-1 is allowed by 
our models presented, if the companion is spinning faster than 
orbital pseudo-synchronization right before the core collapse 
event and the tides exerted on the companion are weaker than 
the nominal theoretical values. 

9.4. Constraints on the BH Kick Direction 

Since the measured eccentricity of Cygnus X-1 is very low 
(see Table 1), and the orbit has not been circularized much 
since the formation of the BH (see Figure 6), CpostSN has to 
be very small too. As discussed in Section 8, we indeed find 
that CpostSN falls within a range of 0.015-0.022. From the 
orbital dynamics at core collapse, the very low CpostSN might 
shed light on the properties of the natal kick imparted to the 
BH. 

During the core collapse event, CpostSN is affected by the 
amount of mass loss, the direction and magnitude of the na- 
tal kick, as well as epreSN- As mass loss tends to increase the 
eccentricity, a natal kick in the right direction and magnitude 
is needed in order to counteract the effect of mass loss and 
result in a very low CpostSN- It turns out that with the obser- 
vationally inferred constraints on Mbh, M2, Vpec.postSN, and 
the extra constraints mentioned in Section 6, the eccentric- 
ities right after symmetric explosions (i.e. no natal kicks) 
are typically ^ 0.1, and overwhelmingly < 0.45. Also, the 
difference between pre- and post-SN eccentricities is mostly 
< 0.15. These link to the requirement of the small Vpec.postSN- 
Larger amounts of mass loss relative to the total mass of the 
pre-SN system not only leads to higher CpostSN, but also larger 
VpecpostSN, and hence they are not allowed. For the same rea- 
son, the required kick velocity for getting the very low post- 
SN eccentricity also needs to be small. In addition, we note 



that under spe cial condition s, CpostSN could be low even though 
fipreSN is high. iHillsl (Il983h showed that if the symmetric su- 
pernova explosion occurs at the proximity of apastron, CpostSN 
could be ~ for a specific range of mass loss. 

As the eccentricity induced by mass loss is low, natal kicks 
do not have to contribute dramatically to make CpostSN fall 
within the observationally required range. Although natal 
kicks have to be constrained in some directions, that con- 
straint is not super narrow. By extracting data from the win- 
ning binaries, we found that 70% have cos(0) < 0; the distri- 
bution in the azimuthal angle for the kicks of the winning 
binaries does not deviate much from the a priori flat distribu- 
tion. As a result, the requirement of a very low CpostSN does 
constrain the natal kick directions, but the constraint is not 
particularly strong. 
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